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ABSTRACT 

Sunyaev-Zeldovich (SZ) effect from a cosmological distribution of clusters carry 
information on the underlying cosmology as well as the cluster gas physics. In order 
to study either cosmology or clusters one needs to break the degeneracies between the 
two. We present a toy model showing how complementary informations from SZ power 
spectrum and the SZ flux counts, both obtained from upcoming SZ cluster surveys, 
can be used to mitigate the strong cosmological influence (especially that of a%) on 
the SZ fluctuations. Once the strong dependence of the cluster SZ power spectrum 
on (Tg is diluted, the cluster power spectrum can be used as a tool in studying cluster 
gas structure and evolution. The method relies on the ability to write the Poisson 
contribution to the SZ power spectrum in terms the observed SZ flux counts. We test 
the toy model by applying the idea to simulations of SZ surveys. 
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1 INTRODUCTION 

In the last decade, X-Ray observations of galaxy clusters 
have been used extensively to determine the cosmological 
matter density parameter and the amplitude of density fluc- 
tuations (Henry 1997, Bahcall & Fan 1998, Viana & Liddle 
1999). In near future, advent of new mm and sub- mm high- 
sensitivity experiments will open a new window for studies of 
galaxy clusters through the SZ effect surveys. These surveys 
(for example SPT*, ACT^, APEX-Szt) would be capable, 
through wide area coverage and high sensitivity, of detect- 
ing thousands of galaxy clusters up to redshift 2 > 1. The 
cosmological possibilities of such a large data sets are enor- 
mous and will allow to carry out independent estimations of 
the cosmological parameters (Diego et al. 2002, Levine et al. 
2002, Weller 2002, Majumdar & Mohr 2003a,2003b, Hu & 
Kravtsov 2003, Lima & Hu 2004) which could be compared 
with those obtained from other observations (CMB, SNIa, 
Ly-a forest, etc). Specifically, the cluster abundance with 
redshift, would provide constraints on the dark energy 
equation of state parameter w = p/ p (Wang & Steinhardt 
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1998). Since SZ effect is redshift independent, to get the red- 
shifts of the clusters detected in these SZ surveys, one would 
need to have foUowup observations in the optical/IR bands. 
In the absence of redshifts, a large yield SZ survey have 
cluster number counts as a function of SZ flux, M{S) would 
be available. This is also sensitive to cosmological parame- 
ters (Barbosa et al. 1996, 1998, Bartlett 2000) but weaker 
in comparison to the redshift abundance. It is interesting 
then to explore alternative methods which exploit the new 
SZ data sets without any need of z information. This is the 
main essence of the paper. 

In the absence of any redshift information one would 
also be able to estimate the SZ power spectrum from the 
temperature maps in addition to the flux counts. The statis- 
tics of SZ maps have also been well studied, both analytically 
(e.g., Cooray 2000, Zhang & Pen 2001, Komatsu & Seljak 
2002) and through numerical simulations (e.g., Refregier et 
al. 2000, Seljak et al. 2001, Springel et al. 2001, da Silva et al. 
2001, Zhang, Pen & Wang 2002). Comparison made between 
the two approaches (Refregier & Tessier 2002, but also see 
Zhang et al. 2004) show that a simple halo model (Cooray 
& Sheth 2002) description of the clusters reproduces results 
from detailed N-body simulations. Some differences remain, 
especially at very high Vs where presumably the effect from 
substructures creep in. All these studies come to similar con- 
clusions about the potential of using the cluster SZ power 
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spectrum to constrain the background cosmology, especially 
(78 (Bond et al. 2002). Apart from cosmology, it has also 
been shown that with some knowledge of cosmological pa- 
rameters, the SZ power spectrum can be used as a power- 
ful probe of cluster gas physics (Majumdar 2001, Holder & 
Carlstrom 2001, Sadeh & Rephaeli 2004a). It is clear from 
all these studies that the possibilities of doing science with 
the SZ power spectrum are enormous. 

It is, however, well known that cosmological studies 
based on cluster data will show degeneracies in the cos- 
mological parameters and cluster scaling relations. One of 
the best known is the degeneracy between as — fim which 
can only be weakened if one can measure the evolution of 
the cluster population as a function of redshift. Moreover, 
imperfect knowledge of cluster structure and evolution can 
weaken cosmological constraints obtained from these sur- 
veys. In the presence of any such evolution of gas physics 
mimicking cosmological evolution of the cluster number den- 
sity, one would have to resort to other options (Majum- 
dar & Mohr 2003b, Lima & Hu 2004) to tighten the con- 
straints. In the absence of any redshift information, when 
dealing with flux counts or power spectrum, one encounters 
similar dogoneracios between cosmology and cluster physics. 
One can try to break these degeneracies by directly adding 
constraints from flux counts to those from power spectrum 
resorting to time intensive joint-flts to the two data sets. 
This has the advantage in that it not only gives us the con- 
straints on the different parameters but also the correlations 
between them. However, it is instructive to look whether at 
all one can disentangle cosmology from cluster physics be- 
fore resorting to more time consuming techniques. It is in 
this 'exploratory' sense that we try to look at the possibility 
of mitigating the cosmological influence by combining flux 
counts to SZ power spectrum such that it becomes possible 
to probe gas physics. The intention is not to estimate pre- 
cise constraints achievable from upcoming SZ surveys but 
to look at the feasibility of studying cluster gas physics with 
these surveys. This approach towards studying cluster prop- 
erties would be worthwhile since the present uncertainties in 
determining cluster properties from targetted observations 
remains large (like upto many tons of pcrcents for any red- 
shift evolution of the cluster scaling relations) . 

In this work, we propose to combine the cluster num- 
ber counts as a function of SZ effect flux, A/"(>S), and the 
power spectrum of the sample, Ct, in such a way such as to 
dilute the dependency of the SZ power spectrum on cosmo- 
logical parameters. Each depends on the background cos- 
mology and gas physics in slightly different way than the 
other (for example, the SZ power spectrum depends on the 
detailed distribution of the gas in a cluster whereas the SZ 
flux only depends on the total thermal content of the clus- 
ter). Although both Af{S) and Ce depend on the total fluxes 
of the clusters, the latter also depends on the shape of the 
clusters, the size of cluster being inversely proportional to 
its angular distance. Thus the cosmology-gas physics degen- 
eracies are different for the two quantities. Let us note, once 
again, at this point that both Af{S) and Ce are quantities 
where the implicit redshift dependence is integrated over un- 
like that from Moreover, single SZ survey would give us 
both the quantities and hence no extrapolation is needed 
between observations at different wavelengths. 

The rest of the paper is organized as follows: in §2 we 



lay down the basic formalism and then in §3 we describe the 
toy model of the clusters. The standard SZ power spectrum 
and its dependences on different parameters are shown in §4. 
We present the simulations in §5 and in §6 we show the use 
of hybrid power spectrum as a tool for probing gas physics 
and discuss the systematics. Finally, in §7 we discuss and 
summarize the method and results. 



2 THE HYBRID CLUSTER POWER 
SPECTRUM: FORMALISM 

The total SZ power spectrum from a redshift distribution 
of clusters, biased w.r.t. the underlying dark matter distri- 
bution, can be written as the sum of two terms, the 1-halo 
Poisson term and the 2-halo clustering term. In general, the 
clustering term is dominant for £ < 100 (Komatsu & Ki- 
tayama 1999) where the shape of the power spectrum is 
determined mainly by the distribution of the halos. This re- 
gion is driven mainly by the cosmological model. On the 
other hand, the shape of the individual gas profiles deter- 
mine the higher ^-values (around and beyond the peak of 
the total SZ power spectrum). Thus, this part is more in- 
fluenced by cluster physics. Taken together, there remains 
a large amount of degeneracy between the gas physics and 
cosmology. 

The Poissonian contribution to SZ power spectrum, Cf, 
from galaxy clusters can be described as. 



C 



(1) 



where dV{z)/dz is the comoving volume element, 
dn{M, z) jdM is the cluster mass function (Press & 
Schechter 1974, Sheth & Tormen 1999, Jenkins et al. 2001) 
and pi (M, z) is the power spectrum (multipole decomposi- 
tion) of the SZ effect 2D profile of a cluster with mass M at 
redshift z. 

Contribution to the total SZ power spectrum arising 
from the spatial correlations of the clusters, Cf is given by 
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where 6(M, z) is the time dependent linear bias factor. The 
matter power spectrum, Pm{k,z), is related to the power 
spectrum of cluster correlation function Pc{k, z) through the 
bias, i.e., Pc{k,z) = b'^{M, z)D'^{z)Pm{k, z = 0) where we 
adopt h{M,z) = 1 + {iy^{z) - l)/(5c(z) (Mo & White 1996). 
In the above equation D{z) is the linear growth factor of den- 
sity fiuctuation, Sdz = 0) = 1.68 and iy{z) = 5c{z)/a{M). 

As we have already mentioned, future SZ surveys would 
detect many thousands of galaxy clusters giving us a flux 
limited sample A/'(<S) of clusters. Since cluster detection does 
not care about the spatial correlations of the clusters, hence 
A/" (iS) is simply the Poisson distributed sample of clusters. 
It is easy to see that the SZ contribution from this sample 
would give us the 1-halo SZ power spectrum once we are 
careful in evaluating eq. (1). One may then try to rewrite 
eq. (1) in terms oiAf{S). 

To proceed let us write equation (1) in the following 
form. 



© 0000 RAS, MNRAS 000, 1-?? 



Hybrid Q 3 



c 



dS N{S,z)pi{S,z) 



(3) 



where the term N(S, z) is the number of clusters with flux 
S and at redshift z per redshift and flux intervals. In the 
previous equation, instead of the masses of the cluster, we 
have used the measured SZ effect flux (see next section). 
The term pi (5, z) is the power spectrum of a cluster with 
flux S at redshift z. 



PeiS,z) = —fe{S,z) 



(4) 



The function fe{S, z) gives the shape of the power spectrum 

as a function of £. These terms are model dependent and de- 
pends sensitively on the cluster structure and its evolution. 

Combining equations (3) and (4) we get the expression 
for the hybrid power spectrum, C/^, as 



dz N{S,z)fi{S, 



(5) 



The trick is to express N{S, z) as a function of the ob- 
served Af{S) (observed number of clusters with fluxes in the 
interval [S-AS/2, 5 + AS/2]). This is possible if we use our 
understanding of the mass function, dn/dM, of daxk mat- 
ter halos from theory and simulations. One can use scaling 
relations between the mass and SZ flux to rewrite the mass 
function in terms of the observed SZ fluxed. 

Thus, we have 



N{S,z) = A/'(S) 



nth{S) 



(6) 



where nth{S) is the total (theoretical) number of clusters 
with fluxes in the interval [S - AS/2, S + AS/2]. 



nthiS) 



dz 



dV 
dz 



dS 



AS 



dn 
dS 



(7) 



Equations (5) and (6) condense the main philosophy of 
our method. Build the SZ power spectrum as a function of 
the flux function instead of the mass function (equation 5) 
and use a normalized version of the flux function (equation 
6) such that its integral in redshift is the observed number 
counts, that is: 



Af{S) = / dzN{S, z) 



(9) 



(8) 



Equation (6) is particularly interesting because it con- 
tains the cosmological dependence in both the volume el- 
ement dV /dz and the number counts dn/dS which is con- 
nected to the mass function through the chain rule. 

dn dn dM 
dS " dM~dS 

where dn/dM is the cosmology-dependent mass function. 
The above expression also shows how dn/dS depends on 
the S — M relation. We will see in the next section how this 
relation depends on the T — M scaling relation. Summing 
up, the cosmological dependency of the hybrid power spec- 
trum will appear only in the N{S, z) term but it will be 
diluted because of the normalization by the factor nth(S) in 
equation (6). On the other hand, the dependency on the gas 
physics (or equivalently in the scaling relations) will appear 
in N{S, z) through equation (9), in the S'^ term of equation 
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Figure 1. The solid lines give tlic (normalized) power spectrum 
{fl(M,z)) for various clusters models with (from left to right) 
Tc = 25, 12, 6 and 3 arcmin respectively. Dotted lines are the best 
fitting model given by equations (18) and (19). The bumps in 
the tails are due to the truncated nature of the 2D profile (we 
integrate the SZ distortion in a sphere or R < Ry). 



(5) (because we need to transform masses into fluxes) and 

in the power spectrum of a cluster with flux S at redshift z, 
fi{S, z) as we will see in the next section. It turns out that 
the gas physics dependence is not diluted as much as the 
cosmological one. 

To model the number density of clusters as a function 
of redshift and mass, we use the standard Press-Schechter 
mass function although one is free to choose other mass func- 
tions which render better fits to numerical N-body simula- 
tions. The point to be noted is that in writing the hybrid SZ 
power spectrum we have essentially normalized the cluster 
Poisson power spectrum to the observed flux limited sample 
of clusters. As we will see, this has interesting implication for 
the use of cluster SZ Ci as probes of cosmology and cluster 
physics. The resultant power spectrum has weaker depen- 
dence on cosmology (for example, as) while still retaining 
the sensitivity to gas physics. 



3 MODELING THE GALAXY CLUSTER 

To describe the chistcr structure and evolution, we use a 
simple model which can fit current observations (Diego et 
al. 2001). To describe the SZ profile of a cluster one requires 
only three basic ingredients: the temperature of the clusters, 
their sizes and geometry (i.e radial profile). Based on X- 
Ray observations (Verde et al. 2001, Finoguenov 2001), the 
temperature and size are modeled using the scaling relations 



T = ToM"{l + zf 
Ro = RoM'^{l + zf, 



(10) 
(11) 



where Rv and T are the virial radius and the isothermal 
temperature of the cluster. In the above relations we have 
not only defined cluster structure but have also assumed a 
simple redshift evolution for the cluster structure. The basic 
structure and evolution parameters are taken from Diego et 
al..(2001). 

The third ingredient of the model, the geometry, gives 
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the shape of the 3-D distribution of the electrons in the 
intracluster medium. For this we assume a standard /3-model 
with /3 = 2/3 to describe the spatial distribution of the 
electron density, ne{R), such that 



ne{R) = 



ne(0) 



1 + (R/Rc)^ ■ 



(12) 



This particular model has two free parameters: the core ra- 
dius, Rc, which we take from equation (11) assuming that 
the ratio between the virial and core radius is a constant 
{Rv/Rc = 10). The second free parameter is the normaliza- 
tion. Instead of normalizing with the central electron den- 
sity, ne(0), wo normalize to the total SZ signal (which is an 
observable) from the cluster given by 



TM 



(13) 



where the constant ATo contains the baryon fraction, /t, 
today's mean CMB temperature, Tcmb and the frequency 
dependency of the SZ thermal effect, g{x), such that 



ATo = lQ2g{x)fbTcMB fjK strad" 



(14) 



In the above equations, the temperature is given in 
Kev, the mass in 10^^/i~^Mq, the baryon fraction in 
and the angular diameter distance, Da, in Mpc. The 

/i-dependency of the mass and baryon fraction cancel outs 
with the one in Da- The frequency dependence, g{x) is given 
by 



tanh(x72) 



-4, 



(15) 



where x = i'(GHz)/56.8. 

Although, SZ power spectrum is quoted in //K^, the 
cluster fluxes arc quoted in mjy. Therefore, the correspond- 
ing expression for the total flux would be 



TM 

ST(mJy) = 5„. 



(16) 



As above. So contains all the relevant constants and is given 
by 

So = 4.4 X lOV^/f, mJy (17) 

where f{x) = p(a;)a;'*e"/(e" - 1)^ . 

To calculate the cluster power spectrum (see equation 
3) we also need an expression for the form factor fi{S,z) 
appearing in equation (4). This is obtained by fitting the 
Fourier transforms of SZ profiles of clusters from simula- 
tions whore the mass and redshift of the cluster are changed 
(keeping Rv/Rc constant). A fit of this kind using our model 
to simulate the clusters gives; 



fe{M,z) = 0.9[exp(-^(,.J] + 0.1[exp(-, 
with, 

f _ ,2 1.5/{0.92+10.4r,) 



(18) 



(19) 



where the core radius, Vc, is given in radians. Note that in 
the limit of very largo angular scales (i.e ^ — > 0) we have 
fi{M, z) ^ 1 and the cluster behaves as a point source hav- 
ing a constant power spectrum. In figure (1) we show the 
quality of our fit for various cluster models where we have 
changed the core radius. 
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Table 1. Models used to illustrate the sensitivity of the hy- 
brid power spectrum to the cosmological model and the cluster 
physics. Model A is the fiducial model. The numbers in bold face 
show the parameter which has been changed compared to the 
fiducial model. All numbers are dimensionless except To (Kev) 
and Ro {h~^ Mpc) 



4 COSMOLOGICAL AND PHYSICAL 
DEPENDENCY OF THE STANDARD 
POWER SPECTRUM. 

In order to compute the hybrid power spectrum we need the 
observed curve Af{S). The flux limits of the source counts 
arc determined by the sensitivity limit of the survey giving 
Smin- The maximum flux Smax is determined by the flux 
of the brightest cluster in the catalog. The curve Af{S) is 
usually affected by the selection function of the survey, that 
is, below a certain flux, not all the clusters are detected. 
In our case we will consider that all these points have been 
included and that the final A/'(5) is an unbiased estimate 
of the underlying cluster number counts. To simulate the 
cluster distribution, wo assume that the number of clusters 
in the flux bin i (that is, clusters with fluxes i £ [Si, Si+i]) 
is a Poissonian variable with mean value Si where Si is the 
expected number of clusters in that bin for a given model. 
The error bars assigned to each bin are Poissonian error bars. 
With flux counts determined, we are ready to compute the 
different power spectra. To understand the dependencies of 
Af{S) and Cf /Cf on the underlying cosmological models 
and gas physics, wo look at the following observationally 
motivated models, given in Table 1. 

In the table, model A is our fiducial model. We have 
nine free parameters of which three relates to the underlying 
cosmology (whose changes w.r.t the fiducial model are shown 
in models B,C & D) and the other six parameters define 
cluster structure and evolution (changes to these parameters 
shown in model E,F,G,H,I & J). 

In figure 2 we show two examples of simulated cluster 
number counts for two models (A & C) where we only change 
as- We have also fitted the M{S) from model A (shown by 
smooth dotted hne). One can use the fit to extrapolate the 
model to lower flux limit below the detection sensitivity of 
the surveys. A point to be noted is that the two number 
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Figure 2. Simulated Af{S) curves for an cxpcriiiiciit with the 
characteristics of ACT (150 sq. dcg area, and a few mjy minimum 
flux at 150 GHz) for two models with as = 1-0 (top, model C) 
and (78 = 0.8 (bottom, model A). The other parameters are like 
in model A in Table 1 (except for To for which we use To = 8.5 
keV). The fluctuations at high fluxes follow Poissonian statistics 
as well as the error bars. The dotted line is a flt to model A. 




10 100 1000 10000 



Figure 3. This plot shows the dependence of the standard power 
spectrum given by equations (1) & (2) on cosmological parame- 
ters. The models are described in Table 1 and in the text. 



counts are easily differentiable (for example, they are ~ 8a 
apart at i/ ~ 25mJy and ~ 2cr at f ~ 250ni.]y). Since, the 
flux counts for different cosmologies are easily differentiable, 
it makes sense to normalize the Ct using Af{S). In section 
6, we show how this helps to use Ce to study gas physics. 

Other than the flux counts coming out of the surveys, 
one can use the temperature fluctuations in the survey area 
to construct the total power spectrum, Cf", (given by the 

sum of Cf and Cf)^. The total power spectrum of SZ ef- 
fect from galaxy clusters can be estimated with consider- 
able accuracy using nmltifrequency observations of the CMB 
(Hobson et al. 1998, Bouchet & Gispert 1999, Cooray et al. 

§ To avoid confusion, note that is not the same as which 
we have introduced in the text 



Figure 4. This plot shows the dependence of the standard power 
spectrum given by equations (1) & (2) on cluster structure and 
evolution. The models are described in Table 1 and in the text. 



2000, Maino et al. 2002, Martinez- Gonzalez et al. 2003, De- 
labrouille et al. 2003). 

In figure (3) we show Cf for four models (A, B, C, D) 

which differ in the values of the cosmological parameters. As 
shown numerous times in previous studies, the SZ curves are 
very sensitive to cosmology, especially to as (compare curves 
A, B & C). In fact, the strong dependence of on as has 
been used to argue for high value of crs « 1 based on the CBI 
observations of temperature anisotropics at high multipoles 
(Bond et al. 2002). This value if as is high compared to 
those obtained from primary CMB data, cluster abundance 
or weak lensing studies. The approach of constraining as 
from an estimation of C^, however, has a caveat. It assumes 
that we are sure of the cluster structure and its evolution. 
As we see in next figure, cluster physics can substantially 
infiuence the SZ power spectrum over large range in £- values. 

In figure (4) we look at the influence of cluster gas 
physics on the total SZ power spectrum. In comparison to 
figure (3), it is easy to see that changes in cluster physics 
can be a source of substantial systematics over the entire 
£-range of the SZ power spectrum. Whereas at low £, say 
£ « 1000, changes in gas physics modifies the SZ power 
spectrum slightly less than that due change in cosmology, 
at higher multipoles gas physics has an even stronger effect 
on . Thus, the SZ power spectrum is quite degenerate 
w.r.t changes in cosmology and cluster physics. This degen- 
eracy is strongest between cluster structural evolution and 
cosmology. It is easy to understand the strong degeneracy 
especially at higher ^-values: changing erg strongly affects 
the cluster population and thus changes the amplitude of 
the power spectrum whereas changing cluster physics, es- 
pecially the evolution, affects the high redshift population 
strongly and affects the small angular scales. 

Thus in order to probe either cosmology or gas physics 
from SZ observations one would need to consider the de- 
generacy between the two. As an example, a careful con- 
sideration of the mass-temperature relation (i.e gas physics) 
can yield a weaker dependence of the SZ power spectrum 
on as (i.e ~ cr|) instead of the usual ~ cr| (Sadeh 
& Rephaeli 2004b). The focus of this paper is to look at 
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Figure 5. Simulated 10° X 10° map of SZ survey with ACT Figure 6. Recovered SZ map for constructing A^(5). The clusters 

characteristics at 145 GHz. detected by SEXTRACTOR are marked with ellipses. 



the feasibility of studying cluster physics with SZ surveys. 
To do so one would first need to mitigate the influence of 
as on the SZ power spectrum. As we will show this can be 
done by using the hybrid power spectrum. The first thing 
to realize is that in figure (3), curves A & C (having two 
different erg corresponds to the two curves shown in 

figure (2) whereas, in reality, one would observe only one 
data set giving one specific J\f{S). We will use this fact to 
dilute the strong influence of cts on the hybrid power spec- 
trum. However, before we proceed further, we must show 
that the hybrid power spectrum gives a good description 
of the standard SZ power spectrum. In the next section we 
address this issue. 



5 CLUSTER SURVEY SIMULATIONS AND 
THE HYBRID POWER SPECTRUM 

The central idea of the present approach lies in replacing the 
standard power spectrum with the hybrid power spectrum. 
In the halo model, the total power spectrum can be written 
as the sum of the Poisson contribution and that from the 
spatial correlation of the clusters. Other than the largest an- 
gular scales, the Poisson power spectrum dominates heavily 
over the clustering term (Komatsu & Kitayama 1999). Thus 
above i ~ 200 the effective SZ power spectrum comes totally 
out of contributions from the Poisson distributed objects. 
This total SZ power spectrum would be made out of con- 
tributions from halos distributed over a wide range of mass 
and rcdshifts. At smaller angular scales the halo structures 
have considerable influence. Now, in reality, other than the 
halos there can be diffuse SZ contributions from low density 
filamentary structures. In past simulations (see Springel et 
al. 2001, da Silva et al. 2001, Refregier & Teyssier 2002) it 
has already been shown that the dominant contribution to 
thermal SZ power spectrum comes from halos (or regions 
of high overdensities) such that the analytical halo model 
calculations of SZ power spectrum give reasonable results 
over a wide multipole range. These comparisons were made 
by isolating the halos (or overdensities) in the simulations 
and comparing the power spectrum obtained from these ob- 
jects to the total SZ power spectrum obtained in the simu- 
lations. In actual observations one would only detect halos 
through their observed SZ fiuxes, unlike in the simulation 
where one can identify each halo. Moreover, detecting a SZ 
source would depend on how well we are able to separate 
out the SZ-sky from the observed CMB-sky. The idea of the 
present work is to sec how well we can construct the full 
thermal SZ power spectrum from the flux counts recovered 
in a typical survey. We have already laid down the proce- 
dure to do this in the previous sections and in this section 
we use results from numerical simulations to see whether the 
method works. 

For concreteness, we focus on observations from sim- 
ulated ACT survey which we take to have a flux limit of 
Slim 'a 5 mjy (145 GHz). Depending on this sensitivity, the 



survey would sample only a part of the total cluster popu- 
lation having fluxes S > Sum- We wdll consider an area of 
the sky of 100 sq. deg. which is the expected surveyed area 
after 5 months of observation (1 mouth of integration time). 
ACT is a three band experiment with central frequencies 
at z/ = 145, 225 and 265 GHz with = 20,25,30 GHz 
respectively. The effective resolutions of 1.7, 1.1, and 0.93 
arcmin respectively at these frequencies. The telescope will 
be equipped with a Bolometer array (MBAC) with imaging 
capabilities over an area of 22' x 22' . The planned sensitivity 
of the MBAC per detector are 300, 500 and 700 (/xK sec^^^) 
which for one month of integration time and 100 dcg^ sur- 
veyed area renders sensitivities of a few to around 10 jiK. 

In simulating the SZ observations we include four ingre- 
dients: i) a realization of the primary CMB, ii) full hydro- 
dynamical N-body simulation of the SZ which includes the 
correlation botwceu the chistcrs (from White 2003), iii) the 
expected noise level of ACT in its three bands (145, 225 and 
265 GHz) and iv) the effect of the antenna beam. We neglect 
point sources in the present work and assume all the bright 
point sources to be subtracted out. Note, that presence of 
point sources may be a source of considerable noise for fu- 
ture SZ surveys (White & Majumdar 2003, Knox, Holder & 
Church 2003). 

The CMB map is created from a standard WMAP-likc 
CMB power spectrum (see primary CMB power spectrum in 
figure (7)). The SZ map was kindly provided by M. White 
(from his 'Data archive') and consist of an N-body simu- 
lation normalized to the current estimates of the excess in 
power of the CMB at small scales (White 2003). The angu- 
lar size of the simulation is 10° x 10° with a pixel size of 35 
arcsec. 

After adding the CMB and the SZ maps at the three 
ACT frequencies we include the beam smearing by filtering 
with Gaussian beams of 1.7, 1.1 and 0.93 arcmin FWHM. In 
the last step we incorporate the instrumental noise Eissuming 
it is white Gaussian with an RMS of 6, 10 and 15 //K per 
pixel with a pixel of 35 arcsec. The finally sky map, as would 
observed by ACT, is shown in figure (5). 

At this step there are a variety of methods which can be 
applied to recover the SZ map from these three simulations. 
As an illustration we will use the most naive approach which 
just takes the combination of M145 — M225 as an estimator of 
the SZ effect at 145 GHz. In the previous estimator, M145 is 
the ACT map from the 145 GHz channel and M225 is the cor- 
responding map at 225 GHz but degraded to the resolution 
of the 145 GHz map. Since the frequency dependence of the 
CMB is constant in AT/T units, the CMB disappears com- 
pletely after subtracting M225 to M145. The same happens 
to the kinetic SZ effect which has no frequency dependence. 
This particularly simple estimator gives satisfactory results 
and has the advantage that it does not rely on any assump- 
tion (compared to other component separation algorithms). 
We should note, however, that more elaborated methods 
should render better results since, for instance, we do not 
use the map at 265 GHz which could be useful to increase 
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Figure 7. Shown as a black solid lines is the CMB power spec- 
trum used in the simulation. Red solid line is the total (true) SZ 
power spectrum of the original SZ simulation with no antenna 
convolution. The hybrid power spectrum constructed out of a fit- 
ting to the observed Af{S) is shown with dashed (blue) line. The 
fiducial model is used. The two (green) dotted lines show the 
power spectrum of the recovered SZ map (M145 — M225) (bottom 
green dotted line) and the same power spectrum deconvolved by 
the antenna (upper green dotted line). The antenna- convolved 
power shows an increasing tail at large £'s due to the instrumen- 
tal noise. Note how the hybrid power spectrum is not affected by 
the antenna. 



the detection rate. The final SZ map will show a small bias 
toward higher (more negative) values of the SZ effect due to 
the fact that the 225 GHz map still contains a small positive 
contribution from the thermal SZ. The recovered SZ map is 
a noisy estimate of the true one and we can estimate the 
power spectrum of this noisy SZ map. The result is shown 
as the two dotted lines in figure (7) . The bottom dotted 
line is the recovered SZ power spectrum and the top dot- 
ted line shows the beam-unconvolved power spectrum. The 
true SZ power spectrum (no antenna convolution) from the 
hydro-simulation is shown as a thin solid line. The recovered 
SZ power spectrum (bottom dotted line) can not recover all 
the power at small scales {£ > 3000). The reason is because 
our component separation algorithm (M145 — M225) does not 
include any beam deconvolution and the recovered SZ map 
shows the damping in the small scales due to the antenna 
beam. The increase in power a.t £ > 8000 is due to the in- 
strumental noise which is not affected by the beam smearing. 
One can get a better estimate of the recovered power spec- 
trum combining the 3 frequency maps in an optimal way 
(see e.g Martinez-Gonzalez et al. 2003). 

To construct the hybrid power spectrum we need to get 
an estimate of the cluster flux counts from the noisy recov- 
ered SZ map. We detect the clusters (above 3a) by using 
SEXTRACTOR (Bertin & Arnouts 1996) as our cluster de- 
tection algorithm. In figure (8) we show the cluster number 
counts recovered after SEXTRACTING the clusters from the 
recovered SZ map (triangles) . As a comparison, we also show 
the N{S) curve when we use the same package (SEXTRAC- 
TOR) but on the original SZ simulation (diamonds). The 
drop in counts in the faint end is due to the selection function 
of the survey although this selection function is "method" 
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Figure 8. Comparison of cluster flux counts recovered by SEX- 
TRACTOR (shown with triangles) with actual fiux counts 
found in the SZ simulations (shown with diamonds). The 
line shows a fit to the flux counts. 



dependent. An optimized component separation plus clus- 
ter detection algorithm should render better results; so our 
results should be considered as conservative. Finally, once 
we have extracted the flux counts, we can use the counts 
(or a fit to it as shown in figure (8)) to build the hybrid 
power spectrum. The resultant hybrid power spectrum (for 
our fiducial model A in table 1) is shown as the dashed line 
in figure (7). One can also see that for 100 < £ < 3000, 
the hybrid power spectrum constructed out of J^{S) very 
well mimics the total power spectrum from the hydro sim- 
ulations. The small excess in power at £ £ (100,3000) may 
be due to an overestimation in the fit of the real underlying 
N{S) curve in the bright end. 

In figure (7), notice that at ^ > 3000, the hybrid power 
spectrum does not agree with the actual power spectrum 
(for the cluster model considered here). It is at these high 
^-values that the effect of gas physics is dominant. Thus, 
once the low £ part (say £ < 2000) is fixed by the hybrid 
power spectra (which entirely depends on the observed flux 
counts), one can use the observed high £ part of the ob- 
served SZ power spectra to constrain cluster structure and 
evolution. That is, instead of using the fiducial model for 
the cluster parameters, we chose such cluster parameters so 
as to match the hybrid power spectrum to the true power 
spectrum. 



6 USING THE HYBRID POWER SPECTRUM 
TO STUDY CLUSTER PHYSICS 

In this section we come to the central point of the paper 
which is to present the hybrid power spectrum as a tool 
to study cluster physics. We first demonstrate that the hy- 
brid power spectrum is weakly dependent on the cosmolog- 
ical model. From equation (5), we expect the hybrid power 
spectrum to be less sensitive to as- This is evident in figure 
(9) where we have plotted the hybrid power spectrum while 
changing the cosmological parameters only. The diluted in- 
fluence of the cosmological parameters on the hybrid power 
spectrum is striking when one compares the figures (9) & (3). 
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Figure 9. Sensitivity of the hybrid power spectrum to the cos- 
mological modeL The models are described in Table 1 and in the 
text. The hybrid power is build from the observed curve corre- 
sponding to (78 = 0.8 in figure (2). The fiducial model A is shown 
in bold line. 



Figure 10. Sensitivity of the hybrid power spectrum to the clus- 
ter physics. The models arc described in Table 1 and in the text. 
The normalization is the same as in figure (9) . The fiducial model 
A is shown in bold line. 



At low £'s (large scales) the hybrid power spectrum is com- 
pletely insensitive to the cosmological model. This is because 
at large scales, all the clusters behave as unresolved point 
sources and their power spectrum can be determined directly 
from the cluster flux counts (which is fixed by the observa- 
tions). The power spectrum for the different cosmological 
models only differs at small scales where the redshift distri- 
bution of the clusters sizes (through the angular diameter 
distance) and fluxes (c.f equation 16) plays a part in deter- 
mining the shape of the SZ power spectrum. Also, changing 
the cosmology will change the shape of the N{S, z) surface 
in equation (6) but since N{S, z) is normalized to M{S), the 
dependency on the cosmological model is diluted. The end 
result of the exercise is that the SZ power spectrum loses its 
strong dependence on as and Om- For example, comparing 
models A & C, we see that the two power spectra differ by 
a factor of two or less at £ > 1000, whereas if we had not 
used the extra information from flux counts we would get 
the power spectra differing by factors of four or more for a 
larger £-range. In essence, we have mitigated the influence 
of cosmology on the SZ power spectra making it suitable to 
study gas physics. One can still argue that we can use the 
hybrid power spectrum to constrain the cosmological model 
as it in the case of the standard power spectrum. However, 
in comparison to cosmology the hybrid power spectrum is 
now more sensitive to the cluster physics than in the case 
of the standard power spectrum making the hybrid power 
spectrum a valid tool for studies of cluster physics. 

Next, we explore the capabilities of the hybrid power 
spectrum as a probe of the intracluster medium. We illus- 
trate these capabilities in flguro (10) where we show how the 
hybrid power spectrum behaves once we change the cluster 
structural parameters for a survey which gives a particular 
Af{S) (with cosmological parameters taken from the flducial 
model). A comparison of the sensitivity of the hybrid power 
spectrum to changes in cosmology (figure 9) and changes 
in the cluster parameters (figure 10) shows that the power 
spectrum is more sensitive to the cluster structure and evo- 
lution. On the contrary (see figures 3 & 4), for the stan- 



dard total SZ power spectrum, both cosmology and cluster 
physics have comparable influence on the power spectrum. 
Thus at high ^-values, for the standard power spectrum it is 
difficult to disentangle the effects of cosmology and cluster 
physics whereas for the hybrid power spectrum this degen- 
eracy is considerably weakened. As an example, a change of 
1% in as changes the standard power spectrum by « 8% at 
£ « 4000. A similar change (1 %) in the parameter (f> (tem- 
perature redshift evolution parameter) changes the standard 
power spectrum by 0.82 % at i? ~ 4000 whereas for the hy- 
brid power spectrum the changes are 2.5 % (in as) and 0.85 
% (in 0) respectively. A similar change in the two parameters 
at angular resolutions probed by ALMA (say, £ = 30000), 
changes the SZ power from 7.5 % and 2 %(standard power 
spectrum) to 3.8 % and 1.5 % (hybrid power spectrum) for 
as and (f>. Although these numbers are calculated for changes 
of 1 % in as and 0, we should note that the current uncer- 
tainty in parameters like (j) or -ip arc large (upto 100 % ). 
Whereas, theoretically the self-similar model predicts cer- 
tain redshift dependencies (like = 1 for the M-T relation; 
Kaiser 1986), observationally there is still no concurring evi- 
dence (Ettori etal 2003 and references within). Uncertainties 
in the cosmological parameters are, however, much smaller 
(for example, WMAP gives an uncertainty of about 5% in 
£78 and 15 % in Qm; Bennett etal 2003). Thus, the result- 
ing hybrid power spectrum becomes a powerful probe of the 
cluster physics. 

It is easy to understand the sensitivity of the hybrid 
power spectrum to the underlying cluster physics. To give 
an example, lets take the cluster temperature redshift evolu- 
tion parameter <p. Since the high multipoles are dominated 
by distance clusters (see figure 11), if (f> is larger, distant clus- 
ters will appear hotter and consequently the power spectrum 
at small scales will be larger. This is why model A has more 
power at large £'s than model G. If we take model H, here 
the core radius is smaller for all clusters than in model A, 
consequently the clusters will look more compact and this 
will increase the power at small scales. Model J, on the other 
hand, will have larger cluster sizes than model A (specially 
at high redshifts) so its power will be smaller at small scales. 
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Figure 11. Contribution to the i'-multipolc in the SZ power spec- 
trum as a function of the redshift of the sources Three different 
multipoles are shown, £ = 500, 3000, 10000. The higher the multi- 
pole £, the more dependent is the power spectrum with the high-2 
population. 



curve is not correctly estimated in the high flux regime. As 
an immediate example, we see this error in estimation has 
given us hybrid power spectrum with slightly more power 
than the total power spectrum at large scales (see figure 
(7)). One can try to bypass this systematic by using a fit to 
the Af{S) curve to high fiuxes. However, in order to have a 
good fit it is necessary to have a large number of bright clus- 
ters. On the other hand, a poor sensitivity of the instrument 
will reduce the number of faint clusters detected in the sur- 
vey. The faint clusters will contribute to the hybrid power 
spectrum in the high i regime (where the sensitivity of the 
hybrid power spectrum with the cluster physics is more rel- 
evant) . Moreover, the cluster catalog (and hereby the Af{S) 
curve) will be the product of a processing of the raw data 
(eg compoucnt separation) which usually introduce a se- 
lection function in the cluster catalog. Understanding the 
selection function is fundamental specially at the faint end 
of the catalog where the selection function is expected to 
change dramatically. Usually, one should fit the J^{S) curve 
only using those data points for which the selection function 
is well understood. Then one can extrapolate the fit to lower 
(or higher) fluxes and build Cf from the fitting. 



Since evolution in cluster structure differentiates strongly 
the high redshift population w.r.t the nearby clusters, the 
high multipole moments of d are very sensitive to any evo- 
lution: consequently one can strongly constrain any cluster 
structure evolution with the hybrid power spectrum. 

The power of SZ power spectrum as a probe of cluster 
gas physics has already been realized (Majumdar 2001). Es- 
pecially, it is seen that high l-vidues SZ power spectrum 
is very sensitive to non-gravitational processes (Komatsu 
& Kitayama 1999, Holder & Carlstrom 2001). This feature 
becomes more noticeable when the high mass clusters arc 
subtracted from the surveys. However, the cosmology-gas 
physics degeneracy were not taken into account and the 
cosmological parameters were fixed from other observations 
(say, primary CMB). The novelty of this work is that we 
dilute the cosmological dependence of the SZ power spec- 
trum using complementary information from the same sur- 
vey itself. Moreover, we do not have to worry about cluster 
redshifts. 



6.1 Sources of systematics 

The main source of systematic error in the construction of 
the hybrid power spectrum comes from our partial knowl- 
edge of the Af{S) curve. This partial knowledge may be due 
to the lack of enough statistics or a poor sensitivity in the 
instrument. A survey covering a small area in the sky will 
contain few clusters per flux bin. This will produce a noisy 
estimate of the A/'(5) curve. A small survey area will also 
increase the error bars due to cosmic variance. In particular, 
the latter can have a significant effect in the high flux inter- 
val since we expect bright clusters to be less common than 
the faint ones. A lack of bright clusters in the sample (or sim- 
ilarly an underestimated M{S) curve at high fluxes) will lead 
to an underestimated hybrid power spectrum in large scales 
(small i's). Consequently, by subtracting the hybrid power 
spectrum from the total (observed) power spectrum one can 
get a biased estimate of the correlation term if the J^{S) 



7 DISCUSSION AND SUMMARY 

In the previous sections we have demonstrated that cosmol- 
ogy as well as cluster gas structure and evolution sinmltane- 
ously shapes the thermal SZ power spectrum from clusters 
of galaxies. The power spectrum depends strongly on certain 
cosmological parameters (like as) and on any evolution, if 
present, thus leading to a cosmology-gas physics degeneracy. 
One can take priors from external observations, for example 
flx the cosmological model from primary CMB, and use the 
power spectrum to probe gas physics. However, it is possible 
to combine complimentary information from the same sur- 
vey itself to mitigate the influence of cosmology in order to 
study gas physics. 

To do so, we have formulated a different way of look- 
ing at the SZ power spectrum by constructing the hybrid 
power spectrum. Given that any future large yield SZ sur- 
vey would also detect many thousands of clusters, we would 
have a flux counts of the objects detected in the survey. At 
the same time, temperature fluctuations from all the survey 
pixels would be used to construct the total SZ power spec- 
trum. We have shown that one can rewrite the Poisson part 
of the SZ power spectrum using the information available 
from the SZ flux counts. Since the Poisson power spectrum 
dominates, in general, over the clustering power spectrum, 
the resultant hybrid power spectrum represents the total SZ 
power spectrum sufficiently well up to £ ~ 2000. We also 
used results from numerical simulations to test and verify 
our analytical results. Since both the SZ power spectrum 
and the flux counts represent the same underlying cosmol- 
ogy, the hybrid power spectrum is, in essence, normalized to 
the background cosmology. Once the cosmology-gas physics 
degeneracy is diluted, we axe left with modeling of the gas 
physics to understand the observed total SZ power spectrum 
at higher multipoles. 

We have used a simple model of the cluster gas structure 
and evolution as well analytical results for the mass func- 
tion in our analytical model. This is sufficient for the present 
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work which is exploratory in nature. Needless to say, prac- 
tical application of the method would need better modeling 
of the complex cluster structure using N-body hydro simula- 
tions or with more complicated analytical models (suitable 
to fit simulation results). Cosinological simulations of large 
yield of clusters mimicking upcoming surveys are still naive 
in the sense that they still need some sort of phenomeno- 
logicaJ approach to model the gas physics (especially non- 
adibaticity) and has to make simplifying assumptions as to 
the cluster structure. At the same time simulating any non- 
standard cluster evolution is non-trivial and either done too 
simplistically or is put in by-hand in the simulations. It is 
this very cluster structure and evolution that wo arc able to 
probe by studying the SZ power spectrum at high ^-values. 
Note, that in figure (7), at multipolcs (£ > 2000) the hy- 
brid power spectrum defers from the actual power spectrum. 
This, of course, should be the case since the cluster struc- 
ture (c.f. equations (10) & (H)) used in constructing is 
very different from the gas physics used in the SZ N-body 
simulations. However, it must be kept in mind that a part 
of the difference may arise from increase in power at high 
i's due to presence of substructures in the clusters which 
arc naturally captured in any N-body simulations but not 
modeled analytically. 

The ability to construct C/^ depends crucially on our 
capability to get the cluster flux counts. The main source 
of systematic error in the construction of the hybrid power 
spectrum comes from our partial knowledge of the J\f{S) 
curve. This partial knowledge may be due to the lack of 
enough statistics or a poor sensitivity in the instrument. 
It is also important to understand that processing of the 
raw data can introduce a selection function in the cluster 
catalog. Usually, one should fit the Af{S) curve only using 
those data points for which the selection function is well 
understood and then extrapolate the fit to lower (and/or 
higher) fluxes. This extrapolation will have little effect in 
the constructed hybrid power spectrum if the survey is large 
enough such that it contains enough bright clusters and if 
its sensitivity is good enough to recover clusters with fluxes 
as low as a few mjy (say, at 145 GHz). Future SZ cluster 
surveys (like ACT/SPT/APEX-SZ) satisfy both those con- 
ditions such that the hybrid power spectrum can be a useful 
tool in these surveys. 

The main conclusion of this work can be summarized by 
comparing figures (9 & 10) with figures (3 & 4). By using 
the hybrid power spectrum instead of the standard power 
spectrum, one can dilute the dependency of the power spec- 
trum on the cosmological model and concentrate on studies 
of the cluster physics. The precision in the determination of 
the cosmological parameters has increased dramatically in 
recent times. However, much uncertainties remain in stud- 
ies of cluster physics. Better handle on cluster structure and 
evolution becomes even more challenging when one takes 
in to account the cosmology-gas physics degeneracy in any 
such studies. The hybrid power spectrum helps to soften this 
problem by diluting rrmch of this degeneracy. 

The prospect of probing gas physics using SZ power 
spectrum would depend on our ability to measure the SZ 
power spectrum as very high /-values. This sub arc-min scale 
has the distinct advantage of the primary CMB contribution 
being negligible. However one has to be careful in eliminat- 
ing other possible sources of secondary CMB anisotropics 



that may introduce further systematic uncertainties. In mul- 
tifrequency experiments, an appropriate SZ reconstruction 
algorithm should be able to recover the true SZ power spec- 
trum up to the resolution limit of the experiment 

Finally, all future high sensitivity SZ observations have 
to tackle the noise introduced by unresolved point sources. 
The spectral dependence of the thermal SZ effect would be 
helpful in eliminating many of such systematics. Although 
not studied in this work, another interesting application of 
the hybrid power spectrum is to perform consistency checks 
on the excess in power in single frequency CMB experiments. 
By this we mean that future CMB experiments will observe 
an excess in the power spectrum of CMB fluctuations at / > 
2000. This excess will be due in part to galaxy clusters, and 
in part to non-removed point sources (see figure 2 in White 
& Majumdar 2003). Now, at £ ~ 2000 the hybrid power 
spectrum is a decent representation of the total SZ power 
spectrum without worrying much about cluster physics. If 
there is an estimation of the J^{S) curve for the survey, then 
one can estimate SZ power due to clusters aX, i ^ 2000. Any 
excess power can be interpreted as due to unresolved point 
sources. 

In spite the stringent requirements needed to make SZ 
observations from the upcoming surveys probe of cluster 
physics, the rewards reaped would undoubtedly be great. 
In addition to new high resolution targeted observations, 
the statistical study of the SZ sky from cosmological distri- 
bution of clusters may provide the next leap in the study of 
the structure and evolution of galaxy clusters. 
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Appendix A 

The hybrid power spectrum for X-ray surveys 

The computation of the hybrid power spectrum for the 
case of X-ray surveys can be done following the same pro- 
cedure as in the SZ case. We start with a measured cluster 
number counts, Af{S) (where now the fluxes are given in 
ergs s~^ cm~^). The fiuxes can be connected with the mass 
through the equations; 

" (20) 



4^Df 



(21) 



where Lx is the cluster X-ray luminosity and Di is the lumi- 
nosity distance. The constant Lo will contain all the units 
and the conversion factors (band correction, conversion be- 
tween cts/s and ergs/s). Equation (21) will play the same 
role as the T — M relation in the case of a SZ survey. The 
final ingredient wc need to compute the X-ray hybrid power 
spectrum is the function f{l, M, z) which will differ to the 
one described in section 3. For the case of X-ray images, 
the 2D profile of a cluster will look different since is propor- 
tional to the integral along the line of sight of the density 
squared. The fourier transforms of these 2D profiles can be 
fitted rendering; 

f{l, M,z) = \ (exp(-C;,.J + exp(- V^)) (22) 
with, 

= (23) 



where the core radius, Vc, is given in rads (see Diego et al. 
2003). 
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